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INTRODUCTION 

This study describes the simulation of a non- 
lifting, 2-degrees-of- freedom space vehicle whose 
attitude is controlled by positioning jets directed 
by a digital control system. The simulation is 
performed on the EAI PACE® TR-48 General Pur- 
pose Analog Computer used with the EAI DES-30 
Digital Expansion System ... a general -purpose, 
low cost, digital logic package designed for opera- 
tion with desk- top analog computers. 

The necessary addition of digital functions to the 
all analog simulation of flight dynamics in this 
study makes this problem a true hybrid applica- 
tion, and one normally implemented on large scale 
hybrid systems such as the EAI HYDAC® 2000 
HYbrid Digital-Analog Computing System(l). 

The DES-30 System, although capable of operating 
autonomously for use in digital instruction or 
design, was designed primarily as an expansion to 
the TR-48 analog computer to provide basic hybrid 
capabilities to the small computer facility. It is 
also readily combined with other general -purpose 
analog or digital computers to provide extended 
control functions. As such, it makes available to 
the analog computer in this problem the logic 
functions needed to implement this essentially 
hybrid application. 

EQUATIONS 

A non-lifting, 2-degrees-of-freedom space vehicle 
is described approximately by the following equa- 
tions: 
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where cp is roll angle, P roll rate, 9 pitch angle, 
Q pitch rate, ^ moment of inertia about the x axis, 
Iy moment of inertia about the y axis, I xz moment 
of inertia about the xz plane, M x moment acting on 
the x axis, My moment acting on the y axis. The 
vehicle is operating outside the atmosphere. It 
uses peroxide jets for attitude control. The jets 
are controlled by a control system, part of whose 
operation can be described approximately by the 
equations: 

< 4 > *<p = (<Pc " <P> - K X P 

(5) c 0 = (9 C - 6) - K 2 Q 

where and eg are error signals, indicating the 
discrepancy between the pilot commands <Pq and 
6 C and the actual attitude cp, 9. The terms -K^P 
and -K 2 Q are introduced to provide a velocity 
correction to the error term, to substitute for 
natural damping which is absent. That is, the 
error signals are automatically reduced if the 
vehicle is proceeding towards the desired direc- 
tion with a high velocity. 

These error signals are applied to the jet control 
portion of the attitude control system. The opera- 
tion of the jet control can be described approxi- 
mately by the curve in Figure 1. 
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Figure 1: The Operation of the Jet Control System 
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The jets are operated in on-off fashion, with a 
fixed amplitude but variable on-off time. The error 
signals control the ratio of on to off time within a 
basic cycle of 1 second; that is, the jets are "pulse- 
width modulated". When the error is below a cer- 
tain minimum value, e m , the jets are kept off . 
Between e m and €p 9 the on time of the jet in- 
creases from 5± seconds to 6 2 seconds, in pro- 
portion to the error. If the error exceeds €^ , the 
jets are turned on and left on until the error de- 
creases below 

The error is sampled every 1 second ^, following 
which the jets are allowed to correct in accordance 
with Figure 1; then, the remaining error is 
sampled again, and so forth. 

There are 4 pairs of jets as shown in Figure 2. 
Jets #1 and 6 are 
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Figure 2: The Attitude Control Jet Pairs 

fired to obtain pitching down; 1, 3, 5 and 7 to 
obtain rolling counterclockwise; etc. While the 
complete jet selection logic for the 3-degrees-of- 
freedom case is quite complex, for the purposes 
of this study it has been reduced to: 

To pitch down fire 1,6 (Pair A) 

To pitch up fire 2,5 (Pair B) 

To roll clockwise fire 3,7 (Pair C) 

To roll counterclockwise fire 4,8 (Pair D) 

The vehicle and the control system are to be 
simulated. The simulation is to run in realtime 
or faster (in fact, as fast as possible). 

MAXIMA, CRITICAL QUANTITIES, & PARAM- 
ETER RANGES 

-2 rad/sec ^ P ^ 2 rad/sec I x =2000 slug- ft 2 

I xz = 300 slug-ft2 
-2 rad/sec < Q < 2 rad/sec I y = 1800 slug- ft 2 

t/n the actual system, the error is sampled continuously as long as 
it is below % m . Many other deviations are included in this study 
and won't be mentioned, since the principal purpose of the study 
is tutorial. 



-2 rad < <p < 2 rad 

-2 rad < 9 ^ 2 rad iq =5.0 sec 

K 2 =2.0 sec 

-2 rad/sec 2 ^ P ^ 2 rad/sec 2 

-2 rad/sec 2 ^ Q < 2 rad/sec 2 |m x | =1000 ft-lbs 

| M y | = 800 ft-lbs 

where M x My are the moments due to the jet 
forces (F x = 40 lbs, F =40 lbs) and the effective 
moment arms (l x = 25 ft, ly = 20 ft) . 



METHOD OF SOLUTION 

The analog portion is straight forward. 

The pulse-width modulation arid the jet selection 
logic are the portions best implemented on the 
DES-30. The pulse width modulation is generated 
as in Figure 3. 

A saw-tooth voltage of period 1 second is generated. 
Every time the ramp portion of the saw-tooth wave- 
shape is started, the absolute value of the error 
is sampled and held to the end of the basic 1 sec 
period. This sampled value is compared with the 
rising ramp. As long as the ramp is below the 
sampled value of the error, a constant voltage, 
e*, representing "jet on" is held in the on posi- 
tion; when the ramp exceeds the error, the "jet" 
is turned off . 

The jet selection logic is quite simple and is 
shown in Figure 4. Signal e*, as discussed above, 
is ANDed with signals e + , e~ indicating whether 
the error is in the positive or negative direction; 
the AND gate signals are then the commands C + 
C" to the appropriate jet pair. The jets are repre- 
sented by electronic switches which introduce 
moments in the form of Figure 3 into equations 
(1) and (2). 

The sawtooth voltage in Figure 3 can be generated 
by an EMC integrator. 

The deadzone indicated in Figure 1 can be ob- 
tained in a number of ways. The method used here 
is shown in Figure 5. The signal at the output of 
the OR gate is high whenever e > ^m or whenever 
(c + € m ) < 0, that is, e < - e m : these are the 
conditions describing the deadzone. 
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Figure 3: Pulse Width Modulation Scheme 
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Figure 4: Jet Selection Logic 
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Figure 5: Deadzone Simulation Using Electronic Comparators 



SCALING (Normalized Method) 



The scaled equations are: 
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Figure 6 shows the program. 
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Figure 6a: Vehicle and Control System Dynamics 



CHECKOUT (Debugging) 

GENERAL: The analog static test technique breaks 
down when logic elements are included. In particu- 
lar, in situations such as this one where the feed- 
back loop is digital, a straightforward check of the 
highest derivative is not possible. 

Probably the best way to static check a hybrid 
problem is this. First, as much as possible is 
checked by straight analog techniques. Checking 



the A— >D interface is accomplished by setting up 
some conditions on the analog side and observing 
the anticipated results on Flip Flop, Counter and 
GP indicators. Conversely, to check the D-*A in- 
terface, one sets or resets flip flops manually 
and observes the results on the analog. The pure 
logic sections are static tested by holding the 
analog portion and stepping the digital section 
clock by clock. For this purpose, a "checkout 
list" is very useful (see below). This checkout 
list must be prepared with great care to assure 
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Figure 6b: Digital Controller 



the operator that when an expected result does not 
occur (or an unexpected one shows up), he can im- 
mediately assume a patching error or machine 
malfunction. The preparation of the checkout list 
also forces the programmer to review his logic 
program systematically. More often than not, 
logical errors will be discovered (this is due to 
empirical rule: "No program is ever right the 
first time"). 

If the logic section passes the clock-by-clock 
check, a dynamic check is rarely needed. How- 



ever, whenever possible, a dynamic check should 
be performed. For example, the sawtooth gen- 
erator in this program can be checked dynamical- 
ly by observing its output on a scope, since it is 
independent of the rest of the circuit. The jet 
logic can be checked dynamically by introducing 
a sinusoidal waveform instead of and e ( 



Such tests should be considered early enough in 
the programming so that they may be incorporated 
in and easily switched into or out of the main pro- 
gram. One must plan for checkout , not leave it to 
the last minute after the program is already patched. 
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Figure 6c: Sawtooth Wave Form Generator 



TEST PROCEDURE: The following Test Procedure is 
used: 

1. Clear the DES, place the TR-48 in PS and 
set the pots. Set FF 20 A, 20B, 20C and 
20D manually to set pots 10, 11, 12 and 13, 
respectively. Because of the TR-48 elec- 
tronic switches, these pots should be set 
by observing the outputs of All and A23 
(in RS). FF 21A and 21B should be cleared. 

Place the TR-48 on RS and perform the 
analog static check. SI, S2 are left; S3, 
S4, S5 are right . To read the " second out- 
put" of A01 (and A05, also), reset FF 20A 
and set FF 20B, manually (reset FF 20D 
and set FF 20C, manually, for A05). 

2. Check the A-*D interface as follows: with 
switches S-^ and S2 in the right position, 
adjust P25 to obtain +2 volts at the output 
of integrator 19. Clear the DES. Make 
sure LP 1 is pressed. Then give 1 manual 
clock. Check that GPR flip flops 20A, 20B, 
20C, and 20D behave according to the table 
below. Then change P25 so integrator 19 
reads 0.5 volts, clear, give a clock, and 
check the table again. Reverse the refer- 
ence polarity to P25 and repeat. 
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3. Check the DA interface by setting and re- 
setting flip flops 21 A, 21B manually and com- 
paring the outputs of amplifiers 01 and 05 
against the analog static test list. 

4. Check the logic program by performing the 
operations indicated in the following table: 



Set 



FF 20A 20B 20C 20D 21A 21B 



IND1 

1 
0 
0 

1 



1 
1 

0 

1 

Read 



1 
1 
1 
0 



IND2 

1 
0 
1 
0 



IND3 



IND4 



Dynamically check the ramp generator and 
the jet control logic by switching S^ S 2 to the 
right , pressing LP2, and observing the output 
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Figure 7: TR-48 Potentiometer Assignment Sheet 

of integrator 19 and amplifiers 11 and 23 on 
a scope or a strip chart recorder. Set P14 to 
a desired frequency (not shorter than 1 sec or 
.1 sec inrealtimeor for B= .1 respectively). 
Run the DES with the 1 MC clock. 

6. Return Si_ S 2 to the left; release LP1, S5 to 
left, clear the DES, adjust pots to their run 
values, etc., before going to run by pressing 
SLAVE on the TR-48 and LP 2 ontheDES-30 
(1 MC clock) . 
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Figure 8: TR-48 Amplifier Assignment Sheet 
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